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Abstract 



The solar corona and heliosphere are visible via sunlight that is Thomson-scattered 
off of free electrons, yielding a radiance against the celestial sphere. In this second part 
of a three-article series, we discuss linear polarization of this scattered light parallel and 
perpendicular to the plane of scatter in the context of heliopheric imaging. The differ- 
ence between these two radiances, (pB), varies quite differently with scattering angle, 
compared to the sum that would be detected by a nonpolarizing instrument (B). In 
particular, the Thomson surface defined by 90° scattering angle is a local minimum in 
scattering efficiency for B measurements, but a local maximum in scattering efficiency 
for pB measurements. We describe the polarization properties of heliospheric Thomson 
scattered light and their applications, covering basic scattering physics, signal-to-noise 
considerations, measurement of 3-D object location, background subtraction, and mod- 
eled pB instrument response to particular classes of solar feature. We conclude that 
\ pB measurements of heliospheric material are much more localized to the Thomson sur- 

face than are B measurements, that the ratio pB/B of polarized to unpolarized feature 
brightness can be used to track solar wind features in three dimensions for scientific 
and space weather applications; and that, by providing an independent measurement of 
background signal, pB measurements may be used to reduce the effect of background 
^ . radiances including the stably polarized zodiacal light. 

Subject headings: Sun: solar wind, coronal mass ejections (CMEs), solar-terrestrial 
relations — Solar System: interplanetary medium — Methods: analytical, data analysis 



1. Introduction 

Coronagraphs and heliospheric imagers observe sunlight that has been Thomson scattered off 
free electrons in the cor ona and solar wind. This potential was realized with the invention of the 



coronagraph (|Lvotlll939i ) , and solar wind transients such as coro nal mass ejections (CM Es) have been 



observed with ground-based coronagraphs since the 1950s (e.g. iDeMastus et al.lll973l ). These were 
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accom panied by spacecraft coronagraphs in the 1970s (e.g. iKoomen et al.l Il975l ; iMacQueen et al 
19741 ) and the coronagraph legacy continues to this day. 



The physics by which this light i s scattered ar e well established, with the original theory 
predati ng the discovery of th e electron (ISchusteiill879l ). Other important developments include the 
work of Minnaert ( 1930 ) and Billings ( 19661 ). The latter is the publication most commonly referred 
to when discussing Thomson scattering theory with regard to white light observations. The utility 
of this theory to identify physical properties (such as mass) in solar wind transient phenomen a 
(such as CMEs) obs e rved by coronagraphs is well known. E arly works include iGosling et al.l (|1975l ) ; 



Hildner et al.1 (|l975f ): lRust et al.1 (|1979T ): IWebb et all (|1980T ) 



In the last decade, coronagra phs have been acc ompanied by another type of w hite light observer . 
Heliospheric imagers, first SMEI fevles et alihoosh and then the STEREO/Rls favles et alJbood ). 
observe Thomson scattered light at much larger angles (> 20°) from the Sun , and their ability to 



track transients such as CMEs has been demonstrated (e .g. iTappin et al.ll2004l ; iHoward et al.ll2006l ; 



Webb et al.ll2006l ; lHarrison et al.ll2008l ; iDavis et al,ll2009i ). All current heliospheric imagers observe 
unpolarized Thomson scattered sunlight. 

Although heliospheric imagers use the same scattering physics as coronagraphs, the wide view- 
ing angle leads to significantly different geometry and requires different treatment. For example, 
Thomson scattering becomes much simpler in the heliospheric case because the Su n can be treated 
as a near-point so urce, eliminating the need to carry van de Hulst coefficients (jMinnaertl Il930i ; 
van de Hulstl Il95(j ) when performing scattering/radiance calculations. More immediately, corona- 
graphs are often assumed to operate near the sky plane, but the relevant figure for a wide-field 
imager is the "Thomson Surface" defined by the locus of the point of closest approach to the Sun 
of each line of sight from the observer. That locus is the sphere with diam eter passing between 
the S un and the observer (jVourlidas &: Howard! 120061 ) . Paper I of this series (IHoward &: DeForest 
2012l ) covered the applied theory of Thomson scattering to describe the relationship between this 
broad-field geometry, illumination, and scattering efficiency in unpolarized heliospheric imaging, 
and demonstrates that a fortuitous cancellation yields a broad plateau (the "Thomson Plateau") 
of nearly uniform radiance sensitivity to electron density. In the present paper, II of the series, 
we explore the consequences of Thomson scattering theory for polarized light in the heliospheric 
context, and discuss possible scientific applications of the theory. 

Thoms on scattered polarized light observation has been used since the dawn of the coronagraph 
(jLyotl 1 19331 ) a nd has been used for three dimensional a nalysis of CMEs since shortly after their 



discovery (e.g. iPoland Munrol Il976l : ICrifo et al.lll983l ). The Skylab coronagraph, Solwind, C/P 
on board SMM, and LASCO on board SOHO all had polarizing capabilities. Perhaps because 
polarized cor onagraph imagery , requiring photometry, is harder to work with than is unpolarized 



imagery (e.g. 
recent work (Moran et al 



MacQueenl I1993T) . it has not been fully exp loited in the spaceflight context although 



2010 



de Koning fc Pizzdl201ll ) may indicate a renaissance of polarized 



image exploitation in the corona. 
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In Section [2l below, we discuss heliospheric imagi ng with polarized light in the same context 
as we have previously with unpolarized light (Paper I: iHoward fc DeForestll2012l ). including: basic 
theory; importance of the Thomson surface; a summary description of instrument sensitivity, de- 
tectability, and signal-to-noise; and discussion of two key scientific applications of polarized imaging. 
Section [3] moves to simulations of these effects applied to non-infinitesimal features including CMEs, 
and corotating interaction regions (CIRs). In Section [4] we discuss the applied theoretical results 
and their implications for future missions. 



2. Elementary Polarized Theory 



Thomson scattering theory has been covered at great length by many authors. Finding the 
polarized radiance of a plasma illuminated by a distributed object that subtends significant solid 
angle from the scatter point (e.g. a solar coronal feature only a few solar radii above the photosphere) 
requires geometrical integrals over the range of solid angles of the incoming light, and is even more 
complicated in the general case than is the same calculation fo r unpolarized light. The problem 
has been treat e d in this general case by numerous authors (e.g. Ivan de Hulstl Il950l : iBillinga 119661 ; 
DeForestlll995l : IHoward &z Tappinl 120091 ) as part of coronagraphic image intepretation. But in the 
heliospheric case (as we pointed out in Paper I) the Sun may be treated as a small object, greatly 
simplifying the theory. Here, we re-derive the polarized radianceS formulae in this simpler case, and 
outline the consequences of its form for different aspects of heliospheric imaging through polarizers. 



2.1. Polarized Scattering Basics 



The elementary unpolarized th eory of Thoms on scattering begins with the differential scattering 
cross-section for unpolarized light (|jacksonlll962l ): 



da 

duj 



a t (l + cos 2 x) 



(1) 



where da is the differential cross section for a single electron to scatter unpolarized light into the 
solid angle duj, at = ?" 2 /2 is half the square of the classical electron radius, and \ is the angle of 
scatter, with x = f° r forward scatter. Thomson scattering is absorption and re-radiation (in a 
dipole pattern) of incident electromagnetic radiation. The electromagnetic waves can be resolved 
into two polarized intensity or radiance components: I±or B± of light polarized perpendicular to the 
plane of scatter, and or B\\ of light polarized in the plane of scatter. In local observing coordinates 
with the Y axis aligned radially outward from the Sun, the quantity Iu —I± is the familiar "Q" Stokes 



As in Paper I, we eschew the term "brightness" that has been used confusingly in the literature to represent 
either radiance or intensity, even though we continue using the B and pB abbreviations which originally stood for 
"brightness" and "polarized brightness", respectively. 
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Fig. 1. — Observing diagram shows relevant angles for heliospheric imaging in the context of the 
Thomson scattering geometry (from our Paper I). The line of sight with elongation e passes through 
the scattering site, making an angle of \ with the radial from the Sun. Distance along the line of 
sight is denoted s when measured from the Thomson surface and z when measured from the observer. 



parameter. The two components are scattered quite differently. The perpendicular polarization is 
scattered independently of x, because of the circular symmetry of dipole radiation: 

da 



doj 



2a t , 



(2) 



2a t cos x, 



while the in-plane polarization is projected into the new propagation direction, and the factor of 2 
reflects the fact that the perpendicular polarization is only half of the original unpolarized beam 
while we are now considering a fully polarized beam. For the radial component, the electric field is 
simply projected with cosX) so that the overall intensity is projected scaled by cos 2 x : 

dan 
du 

and it should be clear that Equation ([1]) is simply the average of Equations (J2j) and 
nagraphic and heliospheric imagery viewed through a radial/tangential polarizer are thus quite 
different from the same imagery viewed with no polarizing optics at all. It is most convenient to 
resolve observed radiance features i nto an unpolarized radiance B = + Bu and an "excess polar- 
ized radiance" pB = Bj_ — Bn (e.g. iFisher et al.lll98ll ). Because the propagators are linear, there is 
a ll pB scattering cross-section" made from the difference between the two polarized cross sections 
just as the unpolarized cross-section is made from their sums: 

da T , 



(3) 
Coro- 



°t 1 



cos 2 x) 



a t (sin 2 x) ■ 



Following Paper I, the differential radiance is given by 

dP± -dP\\ , . _ , 



d( P B) 



dcodA 



a t ism x 



TTV 



R 



ds, 



(4) 



(5) 
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where r & and Bq are the radius and radiance (surface brightness) of the Sun, ds is distance along 
a hypothetical line of sight, and dA is surface area normal to that line of sight. The only difference 
between the polarized radiance pB and the radiance B is that the (1 + cos 2 x) term from Equa- 
tions (1) and (2) of Paper I becomes the sin 2 \ term in our Equations (Q| and ([5]). Separating out 
the x-dependent portion of d (pB) echoes Equation (6) of Paper I: 

d(pB) = kTs(£)Gp(x)n e (s,e,a)ds, (6) 

where e is elongation angle, s is distance along a line of sight, a is an azimuthal angle around the 
Sun in the celestial sphere, and n e is the numerical electron density across space. The geometrical 
values are summarized in Figure [TJ The functions kxs an d Gp serve to isolate the \ dependence of 
pB. Those last two terms are given by 

k TS (e) = (B Q a t 7rr 2 Q ) (R sine)" 2 (7) 

and 

G P = sin 4 x, (8) 

where kxs includes the ^-independent (and hence s-independent) parts of d (pB), and the geometric 
factor Gp includes the x-dependent parts. Gp includes both illumination and scatter dependence 
on x, and is sharply peaked around x = 90°. This is in contrast to the unpolarized geometric factor 
G, discussed in Paper I, which has an extraordinariliy broad peak around x = 90°, superosculating 
the f(x) = 1 line there. Figure [2] shows the relationship between G and Gp. 

The two curves G and Gp are tangent at x = 90° because at 90° the Thomson scattered 
light is fully polarized - so pB = B at that scattering angle. The second derivative of Gp is 
strongly negative at the peak, so that in pB measurements there is no Thomson plateau, and pB 
measurements are most sensitive to material close to the Thomson surface ( TS). In particular, pB 



i mage s are moderately well localized to the vicinity of the TS described by IVourlidas &: Howard 



(|2006l ). 



Considering the scattering efficiency factors in G and Gp independently of illumination: the 
scattering efficiency for unpolarized radiance B is minimized on the Thomson surface, while the 
"scattering" effiency for excess polarized radiance pB is maximized on the Thomson surface. 

The polarized excess intensity (pi) of a unit volume of material is the polarized radiance 
integrated over its apparent size. Carrying out the same operations as for Equation (9) of Paper I 
yields: 



d (pi) = P B dn 



B & o- t Ttr 2 & \ 
m J 



sin 6 x 



dN e . (9) 



sin 2 (e) sin 2 (e + x). 

Both d (pB) and d (pi) for dV = lm 3 and n e = lm~ 3 are plotted at constant elongation e and 
at constant heliocentric distance r in Figure [3j versus the out-of-sky-plane angle (or "sky angle") £ 
(with the understanding from Figure [1] that £ = e + x ~ 7T /2). Figure [3] is directly comparable to 
Figure 5 of Paper I. The TS is marked with a single small tick across each line-of-sight plot. The 
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pB-vs-£ curve (upper left) echoes the sharpness of the Gp geometric function in Figure [2J Unlike 
the unpolarized radiance B, there is no local minimum oipB at the TS in the constant-r case (upper 
right). 

Also unlike the unpolarized case, pi has a local maximum for a given size of feature for elon- 
gation angles smaller than 30° (green dashed curve at lower left). At wide elongation angles, 
perspective effects overwhelm the Thomson scattering and illumination efects, and there is no local 
maximum. In all cases, perspective effects skew the location of greatest polarized intensity. The 
20° elongation curve, at lower left, for example, shows a peak excess polarized intensity about 10° 
closer to the observer than is the TS. At constant r, there is a local maximum in excess polarized 
intensity across \ a t ah ex h angles. 



2.2. Sensitivity and Signal-to-Noise Ratio 

As discussed in Paper I, detectability of features visible in pB depends not on their radiance 
but on their total intensity (in this case, total polarized intensity) above a noise floor. Because pB 
and pi are compound measurements (they are formed from the difference between two radiances or 
intensities, respectively), the intrinsic photon noise for pi is greater than for a single measurement 
of I± or In. But in the heliospheric case, intrinsic photonic noise in the intensity of the desired 
Thomson scattered feature is small compared to the large background (and its associated photon 
and other noise sources) that must be subtracted. Even outside the brighter regions of the zodiacal 
light and the galaxy, the brightness is dominated by the starfield (e.g. I Jackson et al. Following 



the analysis in Paper I, we note that the background noise against which features are measured is a 
random variable with approximately constant distribution, with an approximately constant number 
of samples per unit solid angle on the celestial sphere. This means that N = Lf2 ' 5 , where N is the 
noise against which a feature is to be compared, £1 is the solid angle subtended by it, and L is an 
instrument-dependent factor. It is thus possible to calculate a priori how the signal-to-noise ratio 




X (degrees) s/d (unitiess) 



Fig. 2. — Comparison of the pB geometric factor Gp to the unpolarized B geometric factor G 
(derived in Paper I) in both scattering angle \ an d normalized distance s/d shows sharply peaked 
behavior. The sharp peak at 90° is due to the combined local maximum of illumination and the 
sin 2 x dependence of the pB scattering term; see Equations (j4]) and ([5]). 





80 60 40 20 -20 -40 -60 -80 80 60 40 20 -20 -40 -60 -80 

Sky angle ^ (degrees) Sky angle J; (degrees) 



Fig. 3. — Thomson scattering effects localize the pB signal to the TS both at constant e and at 
constant r. At some elongations there is a local maximum of feature excess polarized intensity near 
(but not on) the TS (lower left). See text for discussion. 
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(SNR) varies with geometry for any small visual feature, up to a normalization constant. 

Figure @] shows how the ratio varies for different types of feature and different types of an- 
gular comparison. The left two panels show line-of-sight comparison and the right two panels 
show constant-radius comparison, as in Figure HI The top two panels show behavior of a feature 
with constant volume and mass as it propagates out from the Sun, and the bottom two panels 
show the behavior of a self-similarly expanding feature of constant mass as it propagates outward. 
Real solar wind features likely exist between these two cases. The self-similar expansion with 
constant mass is pessimistic far from the Sun as most features expand laterally but not radially 
as they propagate; and many dense features also accumulate mass through snowplow effects (e.g. 



DeForest. Howard McComasll2012l ). Similarly, the constant- volume expansion is somewhat opti- 
mistic far from the Sun as most features expand laterally to occupy approximately constant solid 
angle relative to the Sun. 

Because the signal (and the SNR) contains factors that are dependent on the specific instrument 
(e.g. aperture and photon efficiency) and on the specific feature being observed, we normalized 
the SNR to be near unity, to display how the ratio varies for a particular class of feature. The 
normalization constants for Figure [3] are numerically the same as the normalization constants we 
used in the analogous figure for unpolarized signal (Figure 6 of Paper I), and the figures may therefore 
be compared directly although the maximum value in Figure |4] is not always unity. Polarized SNR 
is dominated by the interplay of the inefficiency of the pB scattering far from the Thomson Surface, 
and the proximity effect in intensity as the feature grows close to the observer. 

The right-hand plots show, unsurprisingly, that Earth-directed features are not clearly visible 
from Earth in pB alone. Features close to Earth are more readily detectable but that is a perspective 
effect: the peaks in SNR at 60°-80° elongation in all the plots in Figure [4] are due to large apparent- 
size effects in infinitesimal features; the SNR of real features will thus roll off when the feature- 
observer distance z shrinks to the same order as the size of the feature itself: at smaller values of z, 
the small-apparent-feature approximation Q cx z~ 2 no longer holds. 



2.3. Comparison with Unpolarized Scattering: 3-D effects 



Combining polarized and unpolarized imagery, with sufficient SNR, enables location of fea- 
tures in 3-D; th i s technique has been explored ext e nsively in coronagraphs i n recent years (e.g. 



Moran fe Davila 



2004 



Dere. Wang fe Howard 



2005 



de Koning fc Pizzol l201ll ) . Here we develop 



the theory in the simpler heliospheric imaging case. The ratio pB/B (and its feature-averaged 
equivalent, pi /I for a whole feature) is just the ratio of the corresponding geometric factors: 



pB 
~~B 



Gp 
G 



sin 2 x 
1 + cos 2 x ' 



(10) 
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SNR variation: const 8, constant vol. 
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Fig. 4. — Signal-to-noise ratio (SNR) variation for polarized intensity, vs. sky angle £ for four cases: 
at left, constant-e; at right, constant-r At top: constant feature volume; at bottom: self-similar 
expansion. In all cases, the SNR drops more rapidly than in the unpolarized case described in 
Paper I. Observer-directed events have SNRs reduced by a factor of order 10 compared to the ideal 
viewing angle. The curves have been normalized with the same coefficients as Figure 6 of Paper I, 
and are therefore directly comaparable to their counterparts there. 
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Fig. 5. — The ratio of pB to B in small features varies with sky angle £, and is plotted for several 
elogation angles (dashed/dotted curves). Overplotted: two example features exit the Sun at £ = 20° 
(green; nearly in the plane of the sky) and £=80° (purple; just missing the observer), passing through 
elongation angles between 20° and 90°. At each elongation, the feature's pB/B ratio corresponds 
to either the real £ (bold vertical lines with solid circles) or to a ghost location £' that varies with 
elongation (horizontally extended lines with open circles), which is dynamically implausible. See 
Figure [6] for a sketch. 

which has solutions: 



The relationship between £ and pB/B is plotted in Figure [5] for two different features exiting the 
Sun at £ = 20° (near the plane of the sky; green in the figure) and £=80° (nearly directly at the 
observer; purple in the figure). The two branches are on opposite sides of the Thomson surface, 
equally displaced from it along the line of sight. In the coronagraphic case, the two branches 
represent a permanent ambiguity, because e is nearly zero and therefore the TS is in the plane of 
the sky. In the heliospheric case, the curvature of the TS breaks the front/back asymmetry. While 
any one heliospheric measurement of pB/B cannot identify which branch is occupied by a given 
feature, time series of observations can: on the occupied branch, the trajectory is approximately 
inertial and radial from the Sun, while on the opposite branch the inferred "ghost trajectory" includes 
large lateral accelerations, as diagrammed in Figure [6l 




(11) 



or, giving the sky angle £ in terms of observed elongation angle e and pB/B: 




(12) 



Location of features in three dimensions requires propagating noise in the photometry through 
the (pi / I)-to-^ conversion process. Neglecting instrumental noise sources, the noise level of the 
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Fig. 6. — Two feature trajectories show how feature £ measurements relate to elongation angle e. 
The green (£ = 20°) and purple (£ = 80°) features propagate most of 1 AU on opposite sides of 
the Thomson surface. The ratio pB/B (or, equivalently, pi /I) in the feature is consistent with two 
values of £ at each point - but the "ghost trajectory" £'(e) is unphysical. See Figure [5] for a plot of 
pB/B for features along each line of sight in this figure. 

(pi /I) measurement may be calculated by simple quadrature combination of the pi and / noise 
sources. Treating the noise in first order, the feature as compact and uniform-radiance, and the two 
measurements (and noise samples) as independent, gives: 



where denotes addition in quadrature. The resulting noise term can then be propagated through 
Equation (|12p to identify error in £ determination for any given noise level in the images. The 
two terms added in quadrature are the SNRs of the pi and / measurements, respectively. This 
treatment glosses over the fact that pi is a compound measurement and hence incurs two samples 
of the background noise, but that makes no difference because: (a) the photon noise in the signal 
itself is considered to be negligible compared to that of the background; and (b) if pi and / are 
assembled from the sum and difference of two intensities obtained through polarizers (or through one 
rotating polarizer), there are still only two samples of the background photon noise to be considered. 

To understand the behavior of actual 3-D interpretation, we work from the perspective of trying 
to determine the position of a particular feature in radial motion, whose total intensity is measured 




(13) 
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with a given SNR in the unpolarized channel. Figure [7] shows calculated position of an ensemble of 
five small features with moderate SNR, using direct measurement of the pi jl intensity ratio. 

Two of the features in Figure [7J correspond to the two features followed by Figures [5] and [6j 
All of the features were treated as self-similarly expanding, with an SNR of 10, 30, and 100 at 0.8 
AU from the Sun, yielding three sets of error bars on the inferred value of £; these error bars are 
represented with the nested semitransparent curves around each £ trace. To find the SNR at each 
value of r, we used the same formulae as for Figured] (and for the I measurement using the formulae 
for Figure 6 of Paper I). For example, following the green trace we observe that a feature leaving 
the Sun at £ = 20° that happens to have a SNR of 10 at elongation e = 45° can be determined 
from a single pi /I measurement to have £ either between 15° and 25°, or between 67° and 73°. The 
same feature with an SNR of 30 or 100 has similar inferred locations but tighter error bars. 

The shapes of the measurement loci in Figure [7J are dominated by the two branches of the 
£ calculation in Equation (|12p . Since the features propagate radially, the correct branch yields £ 
determinations that are constant as the feature propagates outward through different elongations, 
while the ghost branch (Figure [6]) yields rapidly changing £ angles versus elongation e as the feature 
popagates. The intersection of the two branches occurs at the TS, which is marked with a dashed 
black line. The plot thickens there, primarily because the pB/B ratio is independent of £ at the 
TS. The two branches are symmetric about the Thomson surface because of the symmetry of the 
G and Gp functions around x = 90° (Figure [2]). 

All of the features in Figure [7J show greatly increased uncertainty toward the right of the plot, 
reflecting the fact that the feature grows both fainter and farther from the observer as it propagates 
away from the Sun. The £=0 curve, in particular, has error bars that grow to over ±45° around the 
true measurement at around e =60° , as the SNR drops below unity: in that geometry, the feature 
is 2 AU from the Sun and 2.2 AU from the observer. 

The feature in Figure [7] with the least uncertainty is at £ = 80°, headed nearly directly toward 
the observer. For that geometry, perspective effects oppose the illumination effects (though do 
not exactly cancel them) through the mid portion of the trajectory. Hence, the SNR remains 
high throughout the £ = 80° trajectory, and the error bars remain small except for the intrinsic 
uncertainty near the TS. 

The error bar shape is dominated by the computed SNR of the feature: brighter, more compact 
features yield lower error bars, and fainter, more diffuse features yield larger ones. In particular, 
we have used self-similar expansion with constant N e (total electron count) in the feature to scale 
the error bars as the features propagate outward; this approach may be slightly pessimistic on the 
right-hand side of the plot when compared to actual feature behavior, as discussed in Section 12.21 

We conclude that pB/B measurement of individual features' excess radiance over background 
could be used to infer absolute exit angle from the Sun with precision of a few degrees of angle, with 



SNR values that are achievable with current instruments in B alone (as in lDeForest. Howard Tappin 



20 111 ). While any one image yields ambiguous results with two branches, time series of data reveal 
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which of the two branches is the correct one. Near the TS, such measurements are very imprecise 
because the slope of pB/B with respect to the sky angle £ is near zero - but radially propagating 
features necessarily spend time both near and far from the Thomson Surface. Those events of 
primary interest to space weather prediction - those with £ close to 90° - yield the most precise 
measurements of £, making pB/B measurements particularly interesting for this application. We 
discuss pB/B measurements of large scale features briefly in Section [3] of this paper and at greater 
length (and more rigor) in Paper III of this series (IHoward et al.ll2012l ). 



2.4. Polarization effects and absolute background 



In coronagraphs, pB measurements are used to reduce unwanted background light from sources 
outside the coronagraph. In particular, the terrestrial sky is quite bright near the Sun, and pB mea- 
surements help to reduce that contam ination by subtracting two measurements with equal amounts 
of background light (e.g. iLvotl Il933l ) . In terrestrial coronagraphs, the technique is limited by po - 
larization of the sky light at angles surprisingly close to the Sun (jLerov. Muller fe Poulainlll972l ). 
In space based coronagraphs, the main persistent background is a combination of stray light from 
within the instrument itself and sunlight scattered from microscopic dust grains (the F corona 
and Zodiacal light; we use these terms interchangeably). In the near-Sun field out to ~10-Rq in 
the sky, the F corona is nearly unpolarized, so pB measurements of the corona can provide an 
absolutely calibrate d Thomson radiance o f features in the corona. Unpolarized coronagraphs such 
as SOHO /LA SCO ([B rueckn er et al.lll995T ). and wide- field unpolarized heliospheric imagers such as 
STEREO /HI (jEvles et al.ll2003l ). cannot in principle separate the background signal from a steady 
componen t of the Thomson radiance under stud y, and therefore only report a "feature excess" radi- 
ance (e.g. iDeForest. Howard McComasll2012l ). Here we discuss the feasibility of similar absolute 
calibration of heliospheric Thomson scattered radiance from a polarizing heliospheric imager. 

In a spaceborne heliospheric imager, the background is a combination of stray light from within 
the instrument itself, the unpolarized starfield, and the F corona, which ranges from 10~ 14 B Q to 
lO _13 i?0 at elongations of interest to inner heliospheric solar wind studies in the plane of the ecliptic. 
It is polarized by u p to 20%, with a pole-to-ecliptic radiance variation of about 4x at 90° elonga tion 
from the Sun (e.g. iDumont fc Sanchez! Il976l ; iLeinert. Pitz fc Plancklll98ll ; iLeinert et al.l ll998). In 
other words, there is a background pB signal of order a few xlO _14 i?0 at moderate elongations 
from the Sun, compared to the Thomson-scattering radiance of a few xlO _16 i?0 in typical wind 
features. Thus, even pB measurements require modeling of the diffuse background, just as do direct 
B measurements: individual constructed measurements include a broad diffuse background that 
must be subtracted from the data to yield "feature-excess" radiance in individual features. 

Even with the presence of persistent diffuse background, pB affords techniques for indirect or 
ongoing-calibration methods to improve the modeled background, which could permit detection of 
absolute Thomson radiance based on the improved modeling. In particula r, the Helios probes found 
that the Zodiacal light is extraordinarily stable across solar wind condition (IRichter. Leinert fc Planck 
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Fig. 7. — Inferred sky angle £ vs. elongation e for five small features with moderate SNR, using direct 
measurement of the pi /I intensity ratio. The inferred physical (horizontal) and ghost (diagonal) 
trajectories are solid lines colored by feature. Shapes around each line are £ error bars for self- 
similarly expanding features with three SNRs at 0.8 AU: 10 (wide), 30 (medium), and 100 (narrow). 
The green and purple features are also shown in Figures [5] and [6) The Thomson surface (£ = e) 
is marked as a dashed black line, and curves of constant radius from the Sun are marked as solid 
black curves. 
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1982al lbh and even across the whole solar cycle. iLeinert fc Pita (|1989l ) conclude that the secular vari- 



ation of the Zodiacal light radiance is at most of order 1% (limited by the longevity of the Helios 
photocells) with similarly stable polarization properties. This is useful because, even with as much 
as 1% long-term variation of unpolarized radiance in the Zodiacal light, the stable polarization 
properties can be used to reduce the variation by another 1-2 orders of magnitude. Hence, once 
well characterized with a particular instrument, a Zodiacal light model may be used to produce a 
stable zero-point for pB measurements of Thomson scattering at the few x10~ 17 -Bq level for the 
life of that instrument. 



Other types of wide-field unpolarized background are present in low-Earth orbit (I Jackson et al 



2010), including high altitude orbital speed ram airglow and high altitude aurora (jMizuno et al 



20051 ) . Both effects are thought to be unpolarized as the mechanism for each is direct atomic line 



emission, though to our knowledge neither has been tested. 

We conclude that polarization can remove and/or stabilize at least some important background 
effects, affording long-term remote measurement of absolute density in the solar wind once a suffi- 
cient baseline of measurements is made. While the dominant steady diffuse light source (the Zodiacal 
light) is slightly polarized, it has been measured to be remarkably constant in both intensity and 
degree of polarization; and the dominant variable diffuse light sources in low-Earth orbit (aurora 
and ram airglow) are thought to be unpolarized. 



3. Forward Modeling of Heliospheric Imagery 



In order to gain a ready appreciation of the differences between how solar wind disturbances 
appear in total power and polar ized radiance, we used the generation codes of the Tappin-Howard 
model (jTappin Sz Howard! |2009| ) to compute synthetic sky maps for a number of idealized distur- 
bances. Briefly, a model solar wind density is assumed and the Thomson scattering formulae are 
integrated in a spherical grid. Disturbances are added by multiplying th e relevant cells in the gri d 
by an enhancement factor. These can be bent shells that simulate CMEs ([Howard fe Tappinll2010l ). 
or spiral structures representing corotating interaction regions (CIRs). For the purposes of this 
study, to simplify comparison, we only consider ecliptic plane cuts through those maps. 



3.1. Small-Scale (Blob or Puff) 

The first disturbance considered was a small arc-shaped transient with a thickness of 0.1 AU, 
and a cone angle of 6° (12° total extent) in latitude and longitude, and an enhancment factor of 5 
above the ambient solar wind density. The radiance was computed every 0.1 AU from 0.3-1.2 AU. 
We ran the full simulation for two longitudes in heliocentric observing coordinates, with the center 
of the CME at L=90° (i.e. f = 0°) and L=30° (f = 60°). The locations of the CME are shown 
in Figure [8] (panels a and e). Each colored sector corresponds to a location at which we computed 
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Fig. 8. — Radiance profiles through a small simulated transient. The left hand row (panels a-d) 
shows the case of the transient center at 90° to the Sun-Observer line, while the right row (panels 
e-h) shows the case where it is 30° from the Sun observer line. The top column (a and e) shows 
the locations at which the transient radiance was computed: each color corresponds to a different 
radius; the dashed lines show the central axis of the CME and the Sun-Observer line. The second 
column (b and f) shows the radiance excess as a function of elongation for each radius: the total 
power is solid; polarized power is dashed; the overbars indicate the elongation ranges within which 
the line of sight cuts the Thomson surface inside the transient. The third column (c and g) shows 
the ratio of polarized radiance excess to total radiance excess as a function of elongation. The fourth 
and final column (d and h) shows the ratios of polarized radiance to total radiance. 
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the radiances. In panels b and f we show the excess radiance (i.e. the radiance with the quiet-wind 
values subtracted) for both total power and polarized radiance. In panels c and g we show the 
ratio of the polarized to total radiance excess , and in panels d and h the ratio of polarized to total 
radiance. In all of the plots, the colors of the traces match those in the schematic in panels a and 
e. At elongations where the line of sight intersects the TS inside the transient, an overbar is placed 
on the plots. 

This small transient can be considered to approximate a blob of plasma, and thus provides us 
with a clear view of the differences between observing in total power and observing in polarized 
radiance. 

The most obvious characteristic of the excess radiance plots is the fall-off of radiance as the 
transient moves further from the Sun. It should be noted that although the gradient in elongation 
is much steeper at 90° than at 30°, when we compare similar radii the difference is not particularly 
large, and can be accounted for by the shorter path through the transient at the larger elongations 
for the transient at 90°. 

If we now look at the polarization fraction of the radiance excesses (panels c and g of Figure [8]) 
we see that for the limb transient, the degree of polarization falls off rapidly with elongation. On the 
other hand, for the transient at 30°, the polarization increases to a maximum near 60° elongation, 
and then declines. This can be readily understood if we recall that for scattering at the Thomson 
surface (perpendicular scattering) the scattered light is 100% tangentially polarized. When the 
transient is close to 60° elongation, the intersection of the line of sight and the Thomson surface is 
inside the transient or very close to it, thus the light we receive from the transient is all scattered 
at close to 90°, and so will be highly polarized. For both transient longitudes considered here, the 
signal from the transient at its closest to the Thomson surface has a degree of polarization about 4 
times greater than at its furthest. 

We further see that when the transient is far from the TS, such that its polarization is less than 
that of the quiet solar wind, then the polarization of the total received light is reduced. When it is 
close to the TS the polarization is increased, reaching 80% for the transient directed 30° from the 
Sun-observer line at the time when it crosses the TS. 

3.2. Large-Scale (CME) 

For a large-scale CME, we repeated the simulations for a shell with a longitudinal cone angle 
of 30°. and otherwise identical parameters. The results from this CME are shown in Figure in 
the same format as was used for Figure [8j 

As expected this CME appears much broader in elongation. We do however still see generally 
similar trends with the degree of polarization maximizing near the Thomson surface. It is also 
more evident here than in the case of the small transient that when the leading edge of the CME 




Fig. 9. — Display of the radiances from the large CME simulation. The format is identical to 
Figure [H 
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intersects the TS, the highest degree of polarization (approaching 100%) occurs at the leading edge 
of the CME, thus producing a kind of intrinsic edge enhancement. The trailing edge which is 
remote from the TS, on the other hand, is de-emphasized, in some cases by well over an order of 
magnitude. The CME directed at 30° always extends to the smallest elongations computed (12°) 
when it lies between the Sun and the observer. This is because this CME has a grazing impact with 
the observer. The extension appears as the flat tails in the unpolarized radiance (solid) traces in 
Figure [9F. In the polarized (dashed) traces, on the other hand, the radiance falls off rapidly as the 
scattering approaches pure forward scatter, which is unpolarized. 



3.3. Corotating interaction region. 

The third disturbance we simulated was a spiral structure corresponding to a 400km s -1 CIR. 
This had an extent in longitude of 24°, and increased linearly in density to an enhancement factor 
of 3 at 1.2 AU, beyond which it was of constant density. The results are summarized in Figure [TOl 
in a similar format to Figure [HJ however here we show all elongations from -180° to +180°. The 
computations were carried out at longitudes of the leading edge of the base of 303, 3, 63, 123, 183 
and 243° relative to the observer (the 3° offset is a computational convenience). These correspond 
to approximately -8.5, -4, +0.5, +5, +9.5 and +14 (-13) days from the arrival of the CIR at the 
observer. 

Again a similar trend is seen in that at elongations where the line of sight intersects the TS 
inside the CIR or close to it, the polarization is high (approaching 100%). 

A major difference is that when the CIR is approaching in the East, the polarization at the 
leading edge is low (Figure [TUb and d). This is a direct result of the geometry of a CIR in which 
the tangential contact between the line of sight and the CIR is well beyond the TS whereas for 
most CMEs, when inside 1 AU, the tangent point is generally close to the TS. See Figure 10 of 



Howard &: DeForestl (]2012l ) for plots showing the distance of the CIR leading edge from the TS. 



3.4. Small Feature Intensity vs. Polarized Intensity 

For the small transient, we have also made computations of the total intensity excess (radiance 
excess summed over the whole transient). These were made both for the transient as it propa- 
gates through the heliosphere (constant longitude, varying radius), and also when seen at the same 
elongation but different longitudes. The intensities are plotted in Figure [TO 

When we track the transient at a constant longitude (i.e. as a single transient would be seen 
propagating through the heliosphere: Figure [TT1 a and b) , we see that in all cases, the total intensity 
declines monotonically with radius or elongation. This means that the fall off of radiance is more 
rapid than the increase in solid angle subtended at all radii and longitudes considered. We also see 
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Fig. 10. — Display of the radiances from a simulated CIR. The format is similar to that of Figure [HJ 
with the coloured spirals representing the locations of the CIR at which the simulations were per- 
formed. Since the CIR moves from East (left) to West (right) and is not symmetrical, we show both 
sides of the Sun in the plots. Also the computations are continued all the way to 180° elongation. 
The elongation axis is reversed so that the traces appear in the same orientation as on the sky. Note 
that computations are not performed at elongations less than 12°. 
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Fig. 11. — Plots of the intensity of the small transient as a function of radius at constant longitudes 
(a) , and the same data but plotted against elongation for comparison with Figure [8] (b) ; and as a 
function of longitude at constant elongation (c), for comparison the peak radiance (d). In all panels 
the total power is denoted by solid lines and the polarized radiance by dashed lines, overbars denote 
the approximate region where the line of sight and the Thomson surface intersect inside the CME. 
The panels a and b show the transient at longitudes of 30° (black) , 60° (red) and 90° (blue) . Panels 
c and d show the the transient as seen when it crosses elongations of 30° (black) , 60° (red) and 90° 
(blue). 
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that the polarized intensities fall well below the total intensities away from the Thomson surface. 

The behavior of the intensities of the transient at constant elongation (Figure [Tib ) corresponds 
to how a similar transient would be seen at a particular location in the field of an imager when 
launched at different longitudes. Here we see that for most cases the solid angle subtended by the 
transient is the dominant term, resulting in a monotonic decrease with increasing launch longitude. 
However in the case of the 30° elongation simulation, the polarized intensity does actually increase 
slightly from a longitude of about 20° until the TS. By comparison, the peak radiance (Figure [HH), 
for all cases considered peaks at or near the TS. However for the polarized radiance, the radiance 
is much more sharply peaked, which accounts for the differences of shape in the intensity plots. 



4. Discussion 

We have developed an appropriately-simplified theory of polarized Thomson-scattered imaging 
in the heliosphere, and shown that sensitivity of pB measurements to electron density is well localized 
along the line of sight at the TS, in contrast to the Thomson plateau observed via unpolarized light 
(B). For unpolarized detection, the illumination function and scattering efficiency have equal and 
opposite second derivatives at the TS, leading to the Thomson plateau; in the case of excess polarized 
radiance detection they have equal second derivatives at the TS, leading to a sharper maximum in 
intensity than would be observed via an s— scattering process with no angular dependence. 

The difference in spatial kernel between pB and B enables location of individual solar wind 
features in three dimensions, and we have developed a theory of small feature location including 
the importance of SNR and kinematic effects in determining the precision to which location may 
be measured. The curvature of the TS breaks the front/back asymmetry that hinders attempts to 
accomplish the same thing in coronagraph images near the Sun. Features far from the TS are more 
readily located in three dimensions than are features near the TS, because the pB/B ratio varies 
more with the sky angle £ far from the TS. An instrument that could measure the two relevant 
polarizations with comparable SNR to that of the STEREO /SECCHI imagers could locate the 
exit angle of small features within well under 10°, greatly enhancing interpretation of the remote 
solar wind signal and potentially improving space weather predict ion through direct location of 



Earth-directed features that are hard to measure geometrically (e.g. lLugazll2010l ). 



The pB signal is a differential signal that, by construction, eliminates unpolarized bright fea- 
tures from the field of view; this effect has been used historically in coronagraphs to remove stray 
light and sky light from the Thomson scattering signal. Direct application in the heliosphere is 
more complex because the Zodiacal light is polarized far from the Sun and hence is not removed 
by the pB calculation in the wide-field heliospheric imaging case. However, the Zodiacal light is 
extremely stable both in degree of polarization and in overall intensity, leading to the possibility 
of direct absolute measurements of the Thomson scattering signal using long baselines and in situ 
measurement of the average wind density; the additional stability of the percent-polarization further 



-23- 



reduces the residual background variation below what may be achieved with a single measurement. 
Other diffuse light sources from near Earth include orbital ram airglow and high altitude aurora, 
both of which we anticipate to be unpolarized and hence invisible in pB. 

We have performed simple forward modeling of geometric structures similar to known solar 
wind structures to develop intuition about how a polarized heliospheric imager, if developed, will 
respond to those features. We find that one may expect significant, unambiguous pB signatures of 
the 3-D location of large features, although analysis of such features is more complex than compact 
features whose geometry may be neglected. 3-D location of large scale features such as CMEs is to 
be treated in more detail in the third paper of this series. 

Support for this work is provided in part by the NSF/SHINE Competition, Award 0849916 and 
the NASA Heliophysics program through grant NNX10AC05G. SJT is supported at NSO by the 
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